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ABSTRACT 

It was recently observed in a numerical study on a high order perturbation method 
under heavy fluid loading that a loaded vibrating plate results, not only in the 
classical frequency shift of the in vacuo single resonance (in both the real part 
because of the fluid added mass and the imaginary part because of energy lost by 
radiation), but also in an increase in the number of the resonance. As a result 
of the loading, a single in vacuo resonance of the structure is transformed into a 
multiple resonance. Here we show that this phenomenon is a refinement of the 
Sanchez's classical result where it was established, using asymptotic analysis, that 
in the case of a light loading conditions the scattering frequencies of a fluid loaded 
elastic structure (ie the resonance frequencies) are nearly the real eigenfrequencies 
of the clastic body alone and the complex scattering frequencies of the fluid with 
a rigid solid. A theoretical explanation of the multiple resonances is given using 
classical results on theory of entire functions. It is established that every single 
in vacuo resonance of a simply supported rectangular plate is transformed into an 
infinite number of resonances under fluid-loading condition. 

1 Introduction 

Accurately describe fluid structure interactions increases the computational cost of the 
various numerical methods used. Most numerical methods are leading to the resolution 
of linear systems of simultaneous equations. When the loading conditions are taken 
into account, the size of the matrices involved increases considerably, and they become 
full and frequency dependent. Any method of reducing these drawbacks is therefore 
most welcome; and one of the best methods available is that based on the perturbation 
approach. Asymptotic analysis can be conduced when the loading is light, as in the case 
of a metal structure surrounded by air. This approach consists in introducing a small 
parameter e, which is the ratio between surface mass density and fluid volume density. 
Using perturbation expansions, one can then construct an approximate solution, based 
on the in vacuo eigenmodes, which is not only easier to calculate than the exact solution 
but also leads to a better understanding of the phenomena involved [3] . 

The light Quid loading perturbation method (involving high density structures in con- 
tact with a small density fluid) was recently extended [9J to cases where the perturbation 
parameter becomes large (which corresponds to a light and/or thin structure in contact 
with a high density fluid). During this study multiple resonance processes were observed 
where a single resonance mode can have several resonance frequencies. The spectrum of 
the operator can no longer be described in terms of the discrete spectrum based on a 
countable series of resonance mode/resonance frequency pairs. This typically nonlinear 
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phenomenon is closely linked to the concept of non-linear frequency modes ([2]), where 
the non linearity depends on the time parameters (such as those involved in porous mate- 
rials) rather than on the more usual geometrical parameters. This frequency non-linearity 
is related to problems which are written in the form 0(u)U = S, where u is the unknown 
factor and S is the source term. 0{u) is an operator which depends non linearly on 
the angular frequency u, as in the case of vibroacoustics, where the coupling depends 
non-linearly on the frequency via the Helmholtz equation Green's kernel. 

One shows that this behavior is a general property of at least fluid-loaded simply 
supported rectangular plates. The main results obtained in our previous study [9] are 
recalled in paragraph [2j A theoretical description of multiple resonance in the case of a 
simply supported plate loaded with a non vanishing density fluid is given in paragraph [3j 
It is then established, using classical results on the connection between the order of 
an entire function and the distribution of its zeros, that any resonance mode in this 
structure will have an infinite number of resonance frequencies. Numerical examples of 
a steel plate in contact with air or water are given paragraph SJ The conclusions and 
possible extensions of this work are presented in paragraph 

2 Statement of the problem 

Consider a finite elastic structure occupying a domain S, the behavior of which is 
described by a differential operator A (in the case treated in this communication it is 
a simple bi-Laplacian operator). This structure is loaded with a perfect fluid at rest 
extending to infinity. The equation of the fluid-loaded structure is given by the classical 
integrodifferential equation 

AU(M, uj) - p s huj 2 \U(M, uj) - e J U(M', uj)G(M; M', u)dM^j = F(M) (1) 

where e = pf/p s h is a small parameter in the case of a metal structure (volume mass den- 
sity p s and thickness h) in contact with air (volume mass density pj). G(M; M',u) is the 
Green's function of the Neumann problem in the Helmholtz equation outside the surface 
E. F(M) is the source. To facilitate reading, the weak formulation of the problem is intro- 
duced. One defines the radiation impedance @(U, V, uj) = / B J s U(M')G(M; M', uj)V*{M)dMdM' , 
where U(M) and V(M) are functions defined on the surface E, a(U, V*) stands for the po- 
tential elastic energy of the structure, (U, V*) is the usual inner product and p s huj 2 (U, V*) 
is the kinetic energy. Then the weak formulation reads: find U (M) such that for every 
V(M) one has 

a(U, V*) - A m H ({U, V*) - e/3(U, V*, u)) = (F, V*). (2) 

2.1 Resonance modes and eigenmodes for a fluid-loaded struc- 
ture 

The eigenmodes U m (M, uj) and eigenvalues A m (u) (or equivalently the eigenpulsations 
Ci) m (uS) given by A m (u) = pshuj^iu:)) are the non zero solutions of the homogeneous weak 
formulation (j2J), that is find U m (M,uj) and uj m {uj) such that for every V(M) 

a(U m ,V*)-p s hu 2 m (u)((U m ,V*}-eP(U m ,V\u)^ =0, (3) 
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(3(U,V,uj), that characterizes the energy loss by radiation, depends on the frequency. 
Then the eigenmodes, well suited to deal with harmonic regime, are complex and depend 
on the frequency. 

One introduces the resonance modes, that allow to compute the transient regime in 
a natural way; they are complex and do not depend on the frequency [3] . The resonance 
modes U m (M) and resonances pulsations cu m , given by A m = p s hCj^ are the non zero 
solutions of the following problem: find TJ m (M) and Cb m such that for every V(M) 

a(U m , V*) - p s hC? m (jfi m , V*} - e(3(U mj V*, = 0. (4) 

These two kind of modes, that are identical when the loading is not took into account, 
have a strong link as each resonance pulsation is the solution of the resonance equation 

^m(^m) = ^m; (5) 

and each resonance mode is an eigenmode calculated at the corresponding resonance 
pulsation. 

U m (M) =U m (M,iu m ). (6) 
It is worth noting that in the eigenmodes are easier to compute than the resonance modes. 



2.2 High order perturbation method 

When the loading parameter e is small, the eigenmodes and eigenvalues can be calcu- 
lated using a perturbation method. The following perturbation expansion is introduced 
into the weak formulation (J3j 

U m (M) = uW(M) + eU£\M) + --- + e s U^(M) + --- 
A m = A^ + e A« + --- + e s AW + --- 

where u£\m) and A^ } are the eigenmodes and eigenvalues of the elastic structure in 
vacuo. To facilitate the reading, one poses: can then be 

easily established that one has 

to the order 1 in e:A#(w) = A^(3 mm (uj), 



m \ I m 

(0) 



to the order 2 in e: A£> (u>) = Ajff ( (p mm (u)f + TfflfcxW ^ 



A (0) _ X(0)' 



and so on to the higher orders. Next, as previously shown [9], by grouping the terms and 
by a re-summation technique, one shows that the high order perturbation expansion is 
given by 

A (0) 

AmM * i-ei^y (7) 

which remains valid even for higher values of e. It is worth noting that, compared to the 
classical perturbation expansion of the eigenvalues given by A m (u) ~ A„ (1 + ef3 mm (uj)), 
the high-order approximation (J7J) is particularly interesting as it gives very precise results 
for the same numerical cost. 
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The calculation of the eigenmodes and eigenvalues does not requires any particular 
effort and, for a given u, each eigenpulsation is the solution of the following eigenpulsation 
equation 

.-,(0)2 



~ 2 / \ KMm 



(1 - ef3 mm (uj))' v ' 

Equation (jSJ) shows that the calculation of the various eigenpulsations do not requires 
any computational cost but the evaluation of multiple integrals (5 mm {uj). 



2.3 The resonance equation 

By introducing in equation (jSj) the high order expansion of the eigenpulsation given 
by equation (jSJ), the resonance equation (J5J) is approximated as 

u& (1 - e(3 mm (u m )) « £;£>'. (9) 

This shows that the calculation of the resonance pulsations is completely different 
than that of the eigenpulsations since j3 mm {uo) has a rather complicated dependence on 
the frequency and solving equation ([9]) requires the search for the zeros of non-convex 
complex function which is a challenging task. A numerical example given in the paper [9] 
has shown a curious phenomenon of multiple resonance. That particular example has 
showed that the resonance frequency of a particluar mode is not unique at least two 
resonance frequencies could exist simultaneously for a given geometrical configuration. 
To have a better understanding of this non-usual phenomenon the roots of the resonance 
equation (jUj) are now studied in a more general manner using classical results in the 
framework of the distribution of zero of entire functions [SI E]- in the following the 
resonance equation (jUj) is written as f mm (z): 

r m {z) = z 2 (l-er m (z))-zl m) (10) 
where z = Cj m and z 0m = u>m ■ The zeros of f mm (z) are the resonance pulsations. 



3 Multiple resonances 

3.1 Methodology 

Now all the problem focuses on the existence and the counting of the roots of f mm (z). 
More or less, the theory of the zeros of entire functions (analytic functions in the whole 
complex plane) is based on the Great Picard's Theorem [TTJ which states that if f{z) 
has an essential singularity at a point z w then on any open set containing z w , f(z) takes 
on all possible values, with at most one exception, infinitely often. The exception is 
needed as showed by the function exp(l/z) that has an essential singularity at z = 
but never attains as a value. The second point is that every entire function is either 
a polynomial or has an essential singularity at infinity. Now if one can prove that the 
resonance function f mm (z) is an entire function, that is not a polynomial and that the 
value is not the exception needed by the Picard's theorem, then the resonance function 
f mm (z) has an infinite number of zeros. Then for every mode, an infinite number of 
resonance exists. 



063 - 4 



Now to characterize complex function f(z) as an entire function, the simplest way 
is to study its development as a power series of the form f(z) = ^™I^°c n 2; n . When 
lim^oo ^/c^ = 0, f(z) is an entire function. To show that it is not a polynomial, 
one can study its order. Let us recall that if f(z) is an entire function, the function 
M(r) = max| 2 | =r \f{z)\ will increase indefinitely together with r. The rate of growth 
of M(r) can be characterized by comparing it with function exp(r). By definition, the 
order (of growth) p is the smallest number such that M(r) < expr p+r? for every rj > 
and for r > r (r]). From this, it can be easily seen that the order is given by the 
formula p = limsup r ^ 00 In In M(r)/ In r ; as examples, a polynomial in z is of order 0, 
and the exponential function is of order 1. The order p of an entire function depends 
on the asymptotic behavior of the coefficients of its series expansion. One has 1/p = 
lim n ^ 00 inf(ln(l/|c n |)/(nlnn)) QllJ); or in other words, if c n oc l/{n\) a then p = 1/a. 

One can simplify again the problem under consideration. The resonance function can 
be written as f mm (z) = z 2 — z 2 m — ez 2 j3 mm (z), as a sum of a polynomial [z 2 — z\ m ) of 
order and a function which is the product of a polynomial (ez 2 ) also of order and 
the radiation impedance f3 mm (z) which order is to be determined. The distribution of 
zeros of f mm (z) can be deduced from that of /3 mm (z) as long as e has a non zero value. If 
f3 mm (z) is a polynomial, the resonance function is also a polynomial with a finite number 
of roots. If f3 mm (z) is an entire function (in the case considered below it is of order 1) 
with essential singularity at infinity, then f mm (z) is also an entire function with essential 
singularity at infinity and by the Picard's theorem j mm [^£) can takes on all possible values 
at infinity, with at most one exception, infinitely often. 

The main results of this paper is, as shown below on the particular example of a 
simply supported plate, that the radiation impedance and the resonance function are 
entire functions of exponential type (that is of order one and normal type [5]) that both 
have an infinite number of roots. 



3.2 The resonance function of a simply supported baffled plate 
as an entire function of the frequency 

Now let us consider a thin simply supported rectangular plate (dimensions a x b 
with thickness h) made of isotropic material with Young's modulus E, Poisson's co- 
efficient v and density p s . In a vacuum, the eigenmodes are given by Uml(x,y) = 
2 /\/ab sin mux /a sinmry/b and the eigenpulsations by oOmn = it 2 \/ D / (p s h) [(m/a) 2 + {n/b) 2 ] 
with D = Eh 3 / 12(1— is 2 ). When the plate is baffled, the radiation impedance (3(Umn, Umn, oo) 
is given as 

P{Umn>U^l* ,uj) = / / / / s i n m nx sin niry sin rrnrx' sin niry' — —dxdydx'dy', 

^ Jo Jo Jo Jo d 

(11) 

where d = ^a 2 (x — x') 2 + b 2 (y — y') 2 and c is the speed of waves in the fluid. In the 
following, the pulsation dependence in u is changed into the equivalent wavenumber 
dependence in k = u/c. Since the in vacuo modes involve a double index, the notation 
for the radiation impedance is simplified as ^(UmL, Umn, oS) = f3 mn (k). The resonance 
function is denoted by f mn {k) = k 2 (l — ef3 mn (k)) — k\ mn , with ko m n = Umn/c. 
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It is worth noting that f3 mn (k) and f mn (k) decrease in the upper complex plane (with 
Q(k) > 0, the term exp(tkd) has an exponential decrease when \k\ — *> oo) and conversely 
grow in the lower complex plane, then all the roots of these two functions are located in 
the lower complex half plane. As the time dependency had been chosen as exp(— tut), 
this ensures that the movement of the plate remains finite over time [3]- 

By making a change of variable [TJ, it can be easily shown that (3 mn {k) is given by 

Aab C 1 f 1 e tkD 
r n {k) = / / Q m (x)G n (Y)——dXdY, (12) 

71 Jo Jo D 

where D = \/a 2 X 2 + b 2 Y 2 . The functions G m (x) and G n (y) are given by single analytic 

integrals G m {x) = J sin rmr(x+x') sin mirx'dx' with G m (x) = (smmirx + mTc(l — x) cosmnx) /(2rmr). 

By making a change of variables from rectangular to polar pQ, with < 9 a = arctan(6/ a) < 

7r/2, the radiation impedance (fT2l is written as 

= -- r J i™n(ki 9)d9 - - f' 2 J 2mn (k, 9)d9, (13) 
^ Jo n Je a 

where J lmn (k, 9) and J 2mn (k, 9) are defined, with H mn (R, 9) = G m («) G n (^), as 

a b 

Jimn(k, 9) = H mn (R, 9)e lkR dR, J 2mn (k, 9) = ["*" H mn (R, 9)e lkR dR. (14) 

Jo Jo 

By expanding the exponential as a power series, the two integrals ( fl4l) are given as 

i=oo f^\£ r - g f,=oo f^\(, r - h e 

Jimnjk, 9) = / H mn (R,9)R e dR, J2mn(k,9) = / H mn (R,9)R e dR. 

i=o l - Jo 1=0 l - Jo 

As H mn (R, 9)R l is a continuous integrable function of R in ]0, a/ cos#[ or in ]0, b/ sin9[ 
and of 9 in ]0,7r/2[, the mean value theorem ensures that there exist R± e]0, a/ cos9[, 
R 2 e]O,6/sin0[, 9 1 e]0,9 a { and 9 2 e]9 a ,ir/2{ such that 

f 8a ^ U^y r*/ 2 ^ (ikY 

/ Jimn{k,9)d9 = y j ^ H mn (R 1 ,9 1 )R l , / J 2mn (k,9)d9 =) j ^ H mn (R 2 ,9 2 )R 2 . 

One then obtains f3 mn (k) as an entire function of the variable k: 

r n (k) = -- (lI mn {Ru 9i) £ + H mn (R 2 , 92) £ ■ (15) 

71 V 1=0 t=o ' J 

3.2.1 Order of the radiation impedance 

To calculate the order of the impedance radiation, one studies the comportment of 
the function Mpmn(r) = max| fc | =r \/3 mn (k)\ for r — ► oo. Let us begins with successive 
integrations by parts of the integrals given by equations (fl4")) . Introducing the results 
given in the appendix leads to 

, (h n, -1 , h 2mn (9) h\ mn {9)e^ - h° 3mn (9) d*H mn (R,0) ^ kR 

J im n(k, 0) = ^: + W + j o m W m 

T fir* - 1 ■ W0) h 2 mn (9)e^o -hl mn {9) d*H mn (R,Q) e lkR , ff17 x 

•W*, 0) = - k + — A + | m J^W*) 
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As one studies the comportment of max^^oo \(3 mn (k)\, in the equations (fT6|) and (p3 
the first two power terms in k are negligible compared to the exponential ones when 
| A; | — > oo; similarly, only the first exponential term has a contribution when \k\ — > oo. 
Then 

M pmn (r) =max^ jjh\ mn {e)^d9 + Jjh 2 3mn (e)^a 
Now by expanding the exponentials in equation (fl8l) as power series, one obtains 




M^(r) = max ^ 

' (19) 

As h\ mn (9) and a/ COS6 1 are continuous integrable functions for 6* G [0,6* a ] and ^J mn (#) 
and bj sin# are continuous integrable functions for 9 G [# a ,7r/2]. One applies the mean 
value theorem for the two integrals that appear in equation (fl9l) : there exist two values 
Q\ and 62 such that 

jf = ft -< 9 '> < * < < 20 > 

Let us denote C\ = a/ cos6*i and c 2 = b/ sin 6*2. It is easy to see a < ci < Va 2 + & 2 and 
b < c 2 < Va 2 + 6 2 . Mpm n (r) is then written as 

Mpmn (r) 





1 


max 


|A;| 4 






1 


max 


\k\ 4 





^3mn(^l) 2^ 7j ^ 



£=0 ' e=o 
\h\ mr X0i¥ kCl + hl m M)z lkC2 \ ■ (22) 

The previous result shows that the function Mpmn(r), characterizing the growth of 
the radiation impedance, is the sum of two exponential functions divided by a polynomial 
of the fourth degree. Its order is easy to obtain from classical results on the sum and 
product of entire functions (see for example [IT], chapter VII paragraph 6). There is 
two possibilities for characterizing the growth of Mpmn{r) for r sufficiently large: the 
first is that if a denotes a finite positive number and C a positive constant one has 
Mpmn(r) ~ Cexp(or) which implies that the radiation impedance is of order 1, the 
second is that Mpmn(r) = which implies that the radiation impedance is of order 0. 
From the expressions of the functions h\ mn {9) and h\ mn {9) given in the appendix, it is 
obvious that for a 7^ b and for m 7^ n one cannot have h\ mn {6i)e lkci = — /i| mn (6' 2 )e lfcC2 . If, 
for example c\ > c 2 , then 

ilV,(r) = max^ |^ m .(^i)e lfcci | . (23) 

If a = b and m = n care must be taken; in this case 9 a = tt/4 and a little algebra shows 
that 
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Then the maximum of f3 mm (k) is 

M^ m (r) = max^ |^ mm (^i)e lfcci | • (25) 

Then j3 mn (k) is an entire function of order 1 and the resonance function f mn (k) = 
k 2 (l — ef5 mn (k)) — k 2 mn is an also an entire function of order 1 for all e ^ 0. 

3.2.2 Distribution of zeros of the radiation impedance and of the resonance 
function 

In addition to the order, the relations ( 1221) . ( 1231) and ( 1251) give a more refined infor- 
mation on the distribution of zeros of the radiation impedance. An entire function that 
grows as exp(crr), with < o < oo, is an entire function of order one and normal type 
a. Such functions are known as entire function of exponential type (EFET) [BJ. The 
resonance function and the radiation impedance are EFET. They both possess an addi- 
tional property on the distribution of their zeros. They are functions that belong to the 
Cartwright class C [BJ. A EFET f(z) belongs to the class C if 

- — at < oo, 26 

where the function log + |w| is defined as log + |w| = max(0, log \w\) or equivalently as 
log + \w\ = l/(27r) J^loglw — exp(t8)\d9. A function f(t) bounded on the real axis 
satisfies log + |/(t)| < |/(t)|- To prove that f(z) belongs to the class C it is sufficient 
to prove that \f(t)\ is bounded for every t e 1R and belongs to class C. The radiation 
impedance given by equation (TT3"]) can be majored, for k G M, by: 

\(3 mn {k)\ = £ / (/ H mn {r,e)e tkr dr)d6+ [j" H mn (r,6)e lkr dr } d6 





H mn (r, 9)\dr\de+ \ \H mn (r, 9)\dr)d9) (.27) 



As H mn (r,8) is continuous integrable function of r and 9 the two integrals that appear 
in the right member of equation ([27]) are finite. j3 mn {k) is bounded for every k real and 
belongs to class C. 

Showing that the resonance function f mn (k) = k 2 {l — e(3 mn {k)) + k 2 mn belongs to 
class C requires a more refined estimation as, because of the k 2 term, it is unbounded at 
infinity. Nevertheless, as f3 mn (k) is bounded for every k, the resonance function is also 
bounded for every finite real value k, then for every finite value A, 

+A lo£ + I f mn (t)\ 

8 17 } n dt < oo. (28) 

It remains to estimate j^° lo s jf Wj. ^ As seen in the previous paragraph, the radiation 
impedance for real values of k — > oo has the following asymptotic approximation 




rn{k) ~ W ' ' ^ (% '^^ + ' h lmn{0y k ^ d 6\. (29) 
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The two integrals h{k) = f° a h\ mn {9y k ^d9 and I 2 (k) = jf a h 2 mn {9)e lk ^e d9 that 
appear in this equation can both be approximated for k — > ±00 by the stationary phase 
method [T]: 

Now the radiation impedance equation has the following approximation for k — > ±00 

r-W » ^ ((-l)^^e-f + (-D-^^^J) , (30) 

when — > ±00, /3 mn (k) tends to zero and the resonance function can be approximated 
by f mn {k) ps k 2 . For k — > ±00 one has log + \f mn (k)\ = log \ k 2 \. Then for every value of 
A 

r^ + \ry\ dt w 2 = r-iog^ = 2l± iogA 

(31) 

n_^_QQ jQg-~l~ I f rnrl (t')\ 

Combining results of equations (1281) and (1311) shows that 1+^ < 00. 

Then both impedance radiation and resonance relation belong to the Cartwright class C. 

For < a < it let us denote by n + (r, a) the number of zeros of the function f mn (z) in 
the sector {z : \z\ < r, | arg(z)| < a} and n_(r, a) the number of zeros in {z : \z\ < r, \n — 
arg(z)| < a}. By the theorem of Cartwright and Levinson (theorem 1, page 127 of [6]) 
the density of the set of zeros n±(r, ct)/r has the limit \im r ^ 00 n±(r,a)/r = a/71. This 
result is a very strong one, since not only it ensures that the resonance function possesses 
an infinite number of zeros but also that "almost all" zeros lie in the neighborhood of the 
real axis. 

It is worth noting that the zeros are counted with their multiplicity order. Then the 
zeros of the resonance function appears as conjugate pairs and then must be counted 
twice. Actually, the resonance function f mn (k) = k 2 (l — ef3 mn {k)) + k\ mn can be ap- 
proached for |jfe|^-> 00 by f mn (k) = k 2 (l - e2h\ mm [9 l )/k i e lkci + k 2 mn , then if k is a 
solution of f mn (k) = 0, it is easy to see that —k*, the opposite of its conjugate satisfies 
f mn (-k*) = 0. 



4 Numerical examples 



In this section, some numerical examples of the multiple resonance phenomenon are 
presented. As it is impossible to obtain a closed formula for Ji mn (k,9) and J2mn(k,9) 
given by equations (|14p . a numerically efficient approximation of these functions has been 
constructed. As shown below, this approximation preserves the order properties of the 
involved functions. Since one always has < 9 a < tt/2, the functions to be integrated 
are expanded with respect to the angular variable 9 to the order 1. This gives 

. ia 3 k 3 + 2m 2 7i 2 — iakm 2 ir 2 — (— l) m 2e mk m 2 n 2 , . 

Jlmn(k, 9) ps ji m (k) = a — 1 > , 32 

A(ak — mii) z {ak + rmi) z 

„ „. .,. tb 3 k 3 + 2n 2 n 2 - ibkn 2 ix 2 - (-l) n 2e tbk n 2 n 2 . . 

J2mn{k,9) ^ j 2 n{k) = 6 — . 33 

4(bk — mT) z {nk + nir) z 
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After reporting the two expansions (1321) and (1321) in the radiation impedance given in 
equation (fTBl) . the following approximation of (3 mn (k), denoted by /3 mn (k), is obtained 

P mn (k) = --(jim(k)9 a +j 2n (k)(7r/2 - 6 a )). (34) 

7T 

It is worth noting that the accuracy of the approximation of (3 mn (k) decreases as the 
mode order increases. However, at low orders (say m < 3, n < 3) it remains sufficient. 
It is easy to obtain 

±3t + mn , . 1 - (-l) m . . 

ji m {k = ±rrm a) = a — — , and j lm {k = 0) = a , (35) 

j 2n {k = ±mr/b) = b ±3 [ + W7r , and j 2n (k = 0) = fe 1 - \- . (36) 

Then as ji m {k) and j2n{k) are entire functions of fc, /3 mn (k) is also an entire function of 
A;. It still remains to estimate the order of (3 mn (k). For given m and n and for — > oo, 
one has 

_^_^ m m 2 7r 2 e ?afc — (— l)™72 2 7r 2 e lfefc 
ft^ = 2^ and = 2W ' 

^J mn {k) = 2vr ^-l)™^9 a e iak + (-1)"^ (| - C a ) e lbfc ^) 1 (37) 

It is worth noting that the approximation of the radiation impedance given in equation 
( |34|) or in equation (1371) presents a similar comportment than the exact one: a sum 
of exponential functions divided by a polynomial of the fourth degree; also the lower 
the mode order, the better the approximation (5 mn (k) becomes. To confirm that the 
approximation keeps the order of /3 mn (k) to one, classical results on the sum and product 
of entire functions are used again. First it is worth noting that the entire function 
(—l) m m 2 /a 3 9 a e mk + (—l) n n 2 /b 3 (tt/2 — 9 a )e tbk cannot be identically zero for any value 
of a, b, m or n when \k\ — > oo. To see this more clearly, one can deal separately with the 
cases a = b and a^ b. 

The case a = b. One has 9 a = vr/4 and (3 mn (k) w ^ ((-l) m m 2 + (-l) n n 2 ) e mk ± 
becomes the product of the exponential function e mk of order 1 and the inverse of a 
degree 4 polynomial in k, which is of order 0. /3 mn (k) is therefore an entire function of 
order 1 and of type a. 



then 



The case a ^ b. Let us take for example a > b. Then with \k\ — > oo one has 
max |e mfc | 3> max |e* bfc |; it is worth noting that this is true in the half plane ^s(k) < 0, 
whereas in the half plane Q(k) > 0, $ mn (k) tend rapidly to zero. In the half plane 
Q(k) < 0, one therefore has (3 mn (k) ~ |f ((— l) m m 2 9 a e mk ) p which is also the product 
of the exponential function e mk of order 1 and the inverse of a degree 4 polynomial in 
k, which is of order 0. $ mn (k) is therefore an entire function of order 1 and of type 
a. It is worth noting that if b > a its type is b. In other words, one can assume that 
a = max(a, b). 
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Figure 1: Contour plots of \k 2 (l — ej3 mn {k)) — k 2 mn \ hi the mode m = 2, n = 3. Steel 
plate in contact with air e ~ 0.002 (left) and water e ~ 1.5 (right). 

Both radiation impedance $ mn (k) and resonance relation f mn (k) belong to the Cartwright 
class C and then all, but a finite number of, their zeros must lie in the neighborhood of 
the real axis. 

As shown above the type a of $ mn (k) is given by o = max(a, 6). Then the type 
of f mn (k) is given the type of k 2 (l — ef3 mn (k)) + k 2 mn ; and for e ^ it is given by 
a = max(a, b). By the Cartwright and Levinson theorem, the density of resonances 
lirn r _ +0O (n+ (r, a)/r + n_ (r, a) fr) is then given by 2 max(a, 6) / 7r. To verify this, the sim- 
plest way is to count the zeros of f mn (z) using the argument theorem on the disk C 
with radius K centered on the origin. It is clear from equations fl35l) and fl36l) that 
$ mn (k) and f mn (k) have no pole inside C. The number of its zeros inside this disk, 
in the case of the mode (m — n), is denoted by ni mn {K). It is given by the integral 
n fmn {K) = l/(2nr) § c f mn (k) / f mn (k)dk with k = Ke ld . This integral is computed 
numerically without difficulty. As shown before, the zeros occur in combined pairs of 
negative imaginary part, then each zero must be counted twice; then the density of 
the zeros given by (n(K,a)/K = n + (K,a)/K + n_(K,a)/K has the asymptotic value 
\im K _ 00 n(K,a)/K = 2nj mn (K)/K. 

For example, let us take a steel (E = 200 Gpa, v = 0.3 and p s = 7800 Kg/m 3 ) simply 
supported rectangular plate in contact with air (with density p/ = 1.2 Kg/m 3 and sound 
wave celerity c = 340 m/s) and let the plate dimensions be a = 1 m, b = 0.7 m and 
h = 8.2 cm. It can easily be established that in the mode m — 1, n — 1, one has 
rij n (12) = 2, 71,^(18) = 4 and 72^(30) = 8. In the mode m = 2, n = 3, one has 
rij 23 (19) = 1, nj 23 (20) = 3, "-/ 23 (30) = 7. As the radius K of the disk C becomes large 
one obtains for the density nj mn {K)/K, with e.g. K = 5000 

njii(K) 1 nj 12 (K) _ 1 rip 3 (K) _ 1 n f - 22 (K) 1 
K ~ 3.1407' K ~ 3.1407' K ~ 3.14268' K ~ 3.14268 

This confirms that lim^^ 00 n(i^, a)/K = 2nj mn (K)/K = 2max(a, b)/ir. Similar results 
are obtained if for the fluid water is chosen instead of air or for various plate dimensions. 
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Figure 2: Contour plots of \k 2 (l — ej3 mn (k)) — kl mn \ in the mode m = 2, n = 3. Aluminium 
plate in contact with air e ~ 0.005 (left) and water e ~ 4.5 (right). 

All these numeric values are a validation of the theoretical results obtained in the para- 
graph [3J The resonance relation of the fluid loaded plate has an infinite number of zeros. 
Therefore, while in a vacuum, each mode of a simply supported plate will have only one 
resonance frequency (and its opposite), under loading condition each mode will have an 
infinite number of resonance frequencies. This result is confirmed in figures [1] and [2j In 
these various contour plots, the dark spots that appear correspond to the resonances of 
the simply supported plate, that is the roots of f mn {k) = 0. In figures [U the material of 
the plate is steel with E = 200 Gpa, v = 0.3 and p s = 7800 Kg/m 3 . In figure [2] the plate 
is made of aluminium E = 70 Gpa, v = 0.3 and p s = 2700 Kg/m 3 . The dimensions of 
the two plates under consideration are the same and are given by 1 m x 0.7 m x 8.2 cm. 
While the spectra of both plates are similar in air (excepted in a higher imaginary part 
of all resonances of the aluminium plate) since they both have more or less same ratio 
E/ p s and more or less the same spectrum in vacuo, when loaded by water the difference 
between the two spectra is obvious. It is worth noting that the multiple resonance phe- 
nomenon makes classical modal analysis rather difficult. This could be seen in figure [2] 
where it is easy to se that there is two resonances very close together (around dt(k) ~ 18 
and Q(k) « —3). While in air (or in water for the steel plate) choosing the resonance close 
to the real axis and neglect the resonances with high imaginary part (corresponding to 
rapidly vanishing waves) is natural, in water there is, in the present case, two resonances 
(and then two resonance modes) that have equal importance. 

Now let us return to the initial question of the existence or non-existence of multiple 
resonances in fluid-loaded plates. In every boundary condition, one has to look for the 
resonance pulsations uj mn such that (2)^(1 — e/3™ n (oj)) — Cjmn = 0. Although it has 
not been proved that such a comportment is true under various boundary conditions, it 
seems reasonable to assume that under boundary conditions of all kinds every resonance 
mode in a plate has an infinite number of resonance frequencies as e 7^ (the radiation 
impedance given by relation ffT2]) seems to be an EFET for a large variety of functions 
G m and G n ). 
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5 CONCLUSIONS 



The loading of a rectangular plate with a high density fluid transforms each resonance 
into an infinite number of resonances. These results open up new prospects. In particu- 
lar, it would be interesting to obtain similar estimates on the vibratory modes occurring 
in clamped plates or mechanical structures of other kinds, although these may not be 
easy to obtain. Another question which arises is how to determine under what condi- 
tions moderately coupled systems (such as a steel plate in contact with water) will have 
these properties. Another important aspect worth investigating is the validity of modal 
expansions in the case of resonance modes having two or more resonances frequencies. It 
would also be interesting to develop experimental device which would make it possible 
to test the occurrence of the behavior described here. 
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A Appendix 

The function H mn (r, 9) is given by 



H mn (r,9) 
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Only a few derivatives with respect to r of H mn (r, 9) up to order 3 taken at the points 
r = 0, r = a/ cos 6 and r = b/ sin 9 are not zero. They are given by : 
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